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The study of off-equatorial orbits in razor-thin disks is still in its beginnings. Contrary 
to what was presented in the literature in recent publications, the vertical stability cri¬ 
terion for equatorial circular orbits cannot be based on the vertical epicyclic frequency, 
because of the discontinuity in the gravitational field on the equatorial plane. We present 
a rigorous criterion for the vertical stability of circular orbits in systems composed by 
a razor-thin disk surrounded by a smooth axially symmetric distribution of matter, the 
latter representing additional structures such as thick disk, bulge and (dark matter) halo. 
This criterion is satisfied once the mass surface density of the thin disk is positive. Qual¬ 
itative and quantitative analyses of nearly equatorial orbits are presented. In particular, 
the analysis of nearly equatorial orbits allows us to construct an approximate analytical 
third integral of motion in this region of phase-space, which describes the shape of these 
orbits in the meridional plane. 
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The literature on analytical models for disklike galaxies is vast (see, for instance, 
Ref. [TJ>. Although most disk galaxies contain internal structures such as bars, warps 
and spiral arms,® we can gain a deep insight about the orbital dynamics of stars 
in these systems by modeling them as axisymmetric. This is a standard procedure 
to model, for instance, rotation curves of spiral galaxies. In this framework, a given 
rotation curve is usually fitted by an axisymmetric mass model decomposition into 
a stellar disk, a gaseous disk, a central bulge and a dark matter halo. By assuming a 
constant mass-to-light ratio along the galaxy, it is possible to determine the density 
profiles of the disks by photometry.® 

There is one extra aspect which must be considered in modeling rotation curves 
of disk galaxies: the corresponding mass models must give rise to stable circular 
orbits in the disk’s equatorial plane. If we consider a smooth 3D distribution of 
matter, the stability analysis follows from Taylor expanding the effective poten¬ 
tial up to second order in the vincinity of the circular orbit,® which is stable if the 
epicyclic frequencies k (radial) and v (vertical) are real. However, many galaxy mod¬ 
els in the literature are composed by a razor-thin disk. Although in this case radial 
stability of circular orbits is analyzed as above, vertical stability cannot be analyzed 
by second-order Taylor expansion. This happens because of the discontinuity in the 
vertical component of the gravitational field, due to the surface layer of matter in 
the equatorial plane. In this contribution we briefly present the results of our recent 
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work 3 concerning the vertical stability of circular orbits in Newtonian razor-thin 
disks. After presenting the basics of razor-thin disk modeling, we comment on the 
new vertical stability criterion, which is consistent with the delta-like singularity 
in the density field. Adiabatic approximation permits us to construct an approxi¬ 
mate third integral of motion for nearly equatorial orbits in these systems, which 
describes the shape of the corresponding 3D orbits in the Rz plane. In this way, we 
expect that the formalism presented here can also contribute to the “third integral 
of motion” issue,® as well as to the analysis of orbits in piecewise-smooth dynamical 
systems. Concerning this last topic, we remark that, as far as we know, there is no 
analogue of the KAM theorem which is valid for systems with this kind of discon¬ 
tinuity. However, the literature on disk-crossing orbits presents a large amount of 
numerical evidence in which islands of stability appear in the vicinity of a stable 
fixed point.®® 

The density profile is described as 

p(R,z) = T,(R)6(z) + p s (R,z), (1) 

where is the Dirac delta distribution, p s is smooth and the system has reflection 
symmetry with respect to z = 0. The physical interpretation of the above quantities 
is the following: E(R) is the surface mass density profile of the razor-thin disk and 
p s (R, z) represents all the remaining three-dimensional structures, such as thick 
disk, bulge and (dark matter) halo. For an equatorial circular orbit of radius R 0l 
we found 3 that the condition which must substitute the vertical frequency stability 
condition is 


E(R 0 ) > 0. (2) 

The radial stability condition is the same as in 3D disks, i.e., k 2 (Rq ) > 0. Moreover, 
E > 0 and k 2 > 0 imply (Liapunov) stability of the corresponding circular orbit. 3 
The vertical stability condition presented here is satisfied along the whole region 
of the equatorial plane which contains the surface layer of matter. Therefore, every 
equatorial circular orbit in an axially symmetric razor-thin disk model is vertically 
stable by construction, regardless of the shape of the remaining components of the 
system. 

The procedure considered so far in the literature to analyze vertical stability 
of circular orbits in analytical models of razor-thin disks took into account the 
vertical epicyclic frequency criterion,®^ neglecting the term proportional to S(z). 
This led to the conclusion that some models contain regions with vertically unstable 
circular orbits.®^ 7 A suggested solution to this problem was to consider, in addition 
to the razor-thin disk, thick disks or spheroidal halos in order to guarantee vertical 
stability®However, as described in the preceding paragraph (and discussed in detail 
in Ref. |3|, this is not the case. The addition of 3D distributions of matter does not 
affect vertical stability in such systems. Radial stability, though, is affected by the 
off-equatorial density profile p s . 
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Concerning nearly equatorial orbits close to a stable equatorial circular orbit 
with given z component of angular momentum, adiabatic invariance of the approx¬ 
imate Hamiltonian describing vertical motion near the circular orbit (with R(t) 
acting as an external time dependence in the corresponding potential) gives us the 
following relation 3 for the shape of nearly equatorial orbits in the meridional plane: 

Z(R) _ fZ(R')\ 1/3 

Z(R') V S(i2) J ’ ( J 

where Z(R) is the maximum value of z that the star reaches along the whole tra¬ 
jectory, for each given value of R. This expression represents an approximate third 
integral of motion for the system in the phase-space region under consideration, and 
depends only on the density distribution of the razor-thin disk. The validity of this 
relation was checked numerically in Ref. [3] for a Kuzmin thin disk and a Kuzmin 
disk superposed by a Plummer halo, with good agreement near the equatorial plane 
(see Fig. Q] for orbits in Kuzmin’s potential). This approximate third integral of mo¬ 
tion may be a starting point to apply canonical perturbation theory to integrable 
Hamiltonian systems containing razor-thin disks. 
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Fig. 1. Orbits in Kuzmin’s potential. The solid 
More details can be found in Ref. [3] 
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(red) line corresponds to the prediction of (Hil l. 
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